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Abstract 

In this paper we show that a surface in P 3 parametrized over a 2-dimensional toric variety ST can be represented 
by a matrix of linear syzygies if the base points are finite in number and form locally a complete intersection. This 
constitutes a direct generalization of the corresponding result over P 2 established in [BJ03] and [BC05]. Exploiting 
the sparse structure of the parametrization, we obtain significantly smaller matrices than in the homogeneous case 
and the method becomes applicable to parametrizations for which it previously failed. We also treat the important 
case 2? — P 1 x P 1 in detail and give numerous examples. 
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1. Introduction 



Rational algebraic curves and surfaces can be described in several different ways, the most common being 
parametric and implicit representations. Parametric representations describe the geometric object as the 
image of a rational map, whereas implicit representations describe it as the set of points verifying a certain 
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algebraic condition, e.g. as the zeros of a polynomial equation. Both representations have a wide range of 
applications in Computer Aided Geometric Design (CAGD), and depending on the problem one needs to 
solve, one or the other might be better suited. To give a simple example, the parametric description is better 
for drawing a surface, as it allows to rapidly generate points on the surface, which can then be interpolated, 
whereas an implicit representation is better adapted for testing if a given point lies on the surface, since 
one only needs to check whether the point verifies the algebraic condition that defines the surface. It is 
thus interesting to be able to pass from the parametric representation to the implicit equation. This is a 
classical problem and there are numerous approaches to its solution, see and [CoO"H for good historical 



overviews. However, it t urns ou t that the implicitization problem is computationally difficult. A promising 



alternative suggested in BD07] is to compute a so-called matrix representation instead, which is easier to 



compute but still shares some of the advantages of the implicit equation. Here is the definition. 

Definition 1. Let 34? C P™ be a hypersurface. A matrix M with entries in the polynomial ring K[Tb, . . . , T n ] 
is called a representation matrix of 34? if it is generically of full rank and if the rank of M evaluated in a 
point of P™ drops if and only if the point lies on 34?. 

It follows immediately that a matrix M represents 34? if and only if the greatest common divisor of all its 
minors of maximal size is a power of the homogeneous implicit equation F £ K[Tb, . . . , T n ] of 34?. One major 
ingredient in the construction of such matrices are syzygies. The theory of syzygies has been developed in 
the theoretical context of commutative algebra at the beginning of the 20th century by mathematicians 
such as David Hilbert. However, it was only in the 1990s that the CAGD and geometric modeling commu- 
nity discovered that t he conc e pt of syzy g ies is us eful in their field. Initially unaware of the connections to 
commutative algebra, [SC9fil ]. |SSQK94j . |SGD97^ . and numerous other authors labeled this approach the 



method of "moving curves" (or "moving surfaces") and showed how it can be used to express the implicit 
equation as a determinant. 

In the case of a planar rational curve C given by a parametrization of the form A 1 --■» A 2 , s i— > 
f^fji (s) ) ' wri ere fi G K[s] are coprime polynomials of degree d and K is a field, a linear syzygy (or 
moving line) is a linear relation on the polynomials f±, f%, f^, i.e. a linear form L — h\T\ + /12T2 + h^T^ in 
the variables Ti, . . . , T3 and with polynomial coefficients hi € K[s] such that J2i=i 2 3 hifi = 0- We denote 
by Syz(/) the set of all those linear syzygies forms and for any integer v the graded part Syz(/)j, of syzygies 
of degree at most v. Actually, to be precise, one should homogenize the fi with respect to a new variable 
and consider Syz(/) as a graded module here. It is obvious that Syz(/)j, is a finite-dimensional K-vector 
space and one can easily obtain a basis [Li, . . . ,Lk) by solving a linear system. We define the matrix M„ of 
coefficients of the Li with respect to a K-basis of K[s]„ as 

M v = ( Li L 2 ■ ■ ■ L 

that is, the coefficients of the syzygies Li form the columns of the matrix. Note that the entries of this matrix 
are linear forms in the variables Ti,T^,T^ with coefficients in the field K. Let F denote the homogeneous 
implicit equation of the curve and deg(/) the degree of the parametrization as a rational map. Intuitively, 
deg(/) measures how many times the curve is traced. It is known that for v > d — 1, the matrix M v is a 
representation matrix; more precisely: if v = d— 1, then M v is a square matrix, such that det(M„) = F deg ^\ 
Also, if v > d, then M v is a non-square matrix with more columns than rows, such that the greatest common 
divisor of its minors of maximal size equals F dog ^). 

In other words, one can always represent the curve as a square matrix of linear syzygies. In principle, one 
could now actually calculate the implicit equation. However, it might be advantageous to avoid the costly 
determinant computation and work directly with the matrix instead, as it has the advantage of making 
the well-developed theory and tools of linear algebra applicable to solve geometric problems. For instance, 
testing whether a point P lies on the curve only requires computing the rank of M v evaluated in P. Other 
interesting resul ts using sq uar e matr ix representations directly to solve geometric problems are presented, 
for example, in [ACGS07| or Ma94|, in which intersection problems are treated by means of eigenvalue 
techniques. 
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It is a natural question whether this kind of matrix representation can be generalized to rational surfaces 
defined as the image of a map 



a / . a3 



A — -> A' 

h(s,t) fa(s,t) fa(s,t) 



(s,t) 



fa(s,t)' fa{s,t)' fa(s,t) 



where fi G K[s, t] are coprime polynomials of degree d. In order to put the problem in the context of graded 
modules, one first has to consider an associated projective map 

P i * ( 91 (P) : g 2 (P) : g 3 (P) : 94(P)) 

where is a 2-dimensional projective toric variety (for example P 2 or P 1 x P 1 ) with coordinate ring A 
and the gi G A are homogenized versions of thei r affin e counterparts /j. In other words, S~ is a suitable 
compactification of the affine space (A*) 2 Co03aL Fu93| . In this case, a linear syzygy (or moving plane) of 



the parametrization g is a linear relation on the gi, . . . , g&, i.e. a linear form L = h\T\ + /12J2 + + h{F/± 
in the variables T%, . . . , T4 with hi G K[s, t] such that 

i=l,...,4 

Exactly in the same way as for curves, one can set up the matrix M v of coefficients of the syzygies in a 
certain degree ^, but unlike in the curve case, it is in general not possible to choose a degree v such that M v 
is a square matrix representation of the surface. In recent years, two main approaches have been proposed 
to deal with this problem: 

- One allows the use of quadratic syzygies (or higher-order syzygies) in addition to the linear syzygies in 
order to be able to construct square matrices. 

- One only uses linear syzygies as in the curve case and obtains non-square representation matrices. 

The firs t approach using linear and quadratic syzygies (or moving planes a nd quadr ics) has been treated 
111 |Co03aj for base-point-fr ee homog eneous parametrizations, i.e. 8F = P 2 , and [BCD03I ] does the same in the 
presence of base points. In |AHW 05] . square matrix representations of biho mogen eous parametrizations, i.e. 
ST = P 1 x P 1 , are constructed with linear and quadratic syzygies, whereas [KD06] gives such a construction 
for parametrizations over toric varieties of dimension 2. The methods using quadratic syzygies usually require 
additional conditions on the parametrization and the choice of the quadratic syzygies is often not canonical. 

The second approach, even though it does not produce square matrices, has certain advantages, in par- 
ticular in the sparse setting that we present. In previous publications, this approach with linear syzygies, 
which relies on t he use of the s o-call ed app roximati on com plexes has been developed in the case 2T = P 2 , 
see for example |BJ03| . |BC05{ . and |Ch06l |. and in |BD07l ] for bihomogeneous parametrizations of degree 



(d,d). However, for a given affine parametrization /, these two varieties are not necessarily the best choice 
of a compactification of affine space, since they do not always reflect well the combinatorial structure of the 
polynomials fa, . . . , fa. In this paper we will extend the method to a much larger class of varieties, namely 
toric varieties of dimension 2, and we will see that this generalization allows us to choose a "good" toric 
compactification of (A*) 2 depending on the polynomials fi, ■ ■ ■ , fa, which makes the method applicable in 
cases where it failed over P 2 or P 1 x P 1 and we will also see that it is significantly more efficient and leads 
to much smaller representation matrices. 



The main idea of our method is similar to the one in [BDQ7[. We use a (general) toric embedding to 
consider our domain as a 2-dimensional toric variety contained in a higher-dimensional projective space, 
which we present in Section [2] Contrary to the cited paper, this natural domain will not be in general 
a hypersurface and its coordinate ring will usually not be Gorenstein, which means that we have to give 
new proofs for some of the results in which this property was used. In Section [3] we proceed to establish the 
necessary homological tools and in particular to derive bounds on local cohomology in Theorem llli our main 
technical result. After that, we will see in Section 0] that we can deduce the validity of the approach from 
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previous results to produce an efficient representation matrix for the implicit equation (see Corollary 1 14[) . 
The particular case of bihomogeneous parametrizations of any bidegree is illustrated in Section [5l We then 
show the advantages of our method through several examples in Section [6l After so me co ncluding remarks 



which summarize the scope of the paper in Section[7j an implementation in Macaulay2 M2fl for the important 
special case 2? = P 1 x P 1 is included as an appendix. 

2. Toric embeddings 

Let K be a field. All the varieties considered hereafter are understood to be taken over K. We suppose given 
a rational map 

A 2 --+P 3 
(s,t) ~ (/i : h ■ h ■ h)(s,t) 

where fi € K[s, t] are polynomials. We assume that 

- / is a gcnerically finite map onto its image and hence parametrizes an irreducible surface /cP 3 

- gcd(/i, . . . , f±) = 1, which means that there are only finitely many base points. 

We briefl y intro duce som e basic notion s from t oric geometry. These constructions are investigated in more 
detail in |KD06l . Sect. 2], |Co03bl ]. and |GKZ94l . Ch. 5 & 6]. 



Definition 2. Let p — J2( a p)ez 2 PcP 8 ™^ e K[s,i]. We define the support Supp(p) to be the set of all the 
exponents which appear in p, i.e. 

Supp(p) = {(a, (3) e Z 2 | p a3 £ 0} C Z 2 

The Newton polytope N(/) C K 2 , where / = {fi, fi, fa, fi), is defined as the convex hull of the union 
(J. Supp(/i) in M 2 of the supports of the fi. In other words, N(_f) is the smallest convex lattice polygon in 
R 2 containing all the exponents appearing in one of the fi. Note that our hypothesis that / is generically 
finite implies that N(/) is two-dimensional. Furthermore, let d S N be the biggest integer such that N(/) 
equals 

d-N'(/) = {pi + — +Pd, P,eN'(/)}, 

where N'(/) is a lattice polygon. In other words, N'(/) is the smallest possible homothety of N(/) with 
integer vertices. 

Then N'(/) defines a two-dimensional projective toric variety 8f C P m , as explained in [Co03b| |. where m + 1 
is the cardinality of N' (/) flZ 2 . It is defined as the closed image of the embedding 

(A*) 2 Ap m 

where £ N'(/) n 1? . For example, the triangle between the points (0, 1), (1,0), and (0,0) corresponds 
to P 2 and P 1 x P 1 has a rectangle as polygon. The rational map / factorizes through & in the following way 

(A*) 2 -^p3 (2) 

/ 

9-' 

where g is given by four polynomials g\,...,gi of degree d in m variables. Thus, we have extended the affine 
parametrization / to a parametrization g of over the projective variety 2? 

sr - 9 ^p 3 

P i * ( 9l (P) : . . . : 9i (P)) 
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for which we will adapt the method of approximation complexes. This map induces an application between 
the homogeneous coordinate rings 

K[T u T a ,T 3 ,T4]^A 

Ti i-> Qi(Xo, . . . , X m ) 

where A = K[X , . . . ,X m ]/I(£?) is the homogeneous coordinate ring of which is a domain, as is 
prime. Note that the variables Xk correspond to monomials s z P and the ideal I(^) is the ideal of relations 
between these monomials. The implicit equation of 5? is a generator of the principal ideal ker(/i). We should 
remark that the toric ideals I(^) are very we ll und erstood and there exist highly efficient software systems 



to compute their Grobner bases, for example [4ti2 | 



Instead of N'(/) we could actually have chosen any polygon Q such that a multiple d ■ Q, d € N, contains 
N(.f). I n particular, we could choose N(/) itself, in which case the gi will become linear forms, compare 
KD06I . Sect. 2]. We will see in Section [5] that N'(/) is always a better choice than N(/); for the moment let 



us just state that a smaller polygon leads to a less complicated coordinate ring but to a higher degree of the 
<7i and that the advantages of the former outweigh the inconveniences of the latter. Intuitively, the surface 
3? should be understood to be the smallest compactification of (A*) 2 through which the map / factorizes, 
so in a way it respects the geometry of the map best and is a natural candidate. However, we will see in 
Example [T8l that in some cases there are better choices than the canonical choice N'(/). 

2.1. The combinatorial structure of the ring A 

We can describe the ring A in a more combinatorial way, which will enable us to study its properties in 
more detail. Let C be the cone generated by the polytope N'(/). i.e. the rational cone over N'(/) x 1 C R 3 . 
Then C fl Z 3 equals the union of C n for n £ N, where 

C n = {(i,j,n)\(i,j)e(n-N'(f))nZ 2 }CZ 3 

which means that at each height n we have a homothety of N'(/) by a factor of n. In particular, we identify 
C\ with N'(/) fl Z 2 . Since we are dealing with polygons in dimension two, it holds that 

(n-N'(/))nZ 2 =n-(N'(/)nZ 2 ). 

This property is called normality. We should note that these considerations are no l onger t rue in higher 
dimensions. This is because in dimension > 3 there exist non- normal lattice polytopes MS05L Ex . 12.6]. In 



fact, the study of normality of smooth lattice polytopes is a subject of current research HHM07] 
As an illustration, consider the following picture of the cone C: 




Now we can associate an affine semigroup ring K[C] to this cone: one takes the K-vector space freely 
generated by the elements of C fl Z 3 and equips it with a natural multiplication, which is induced by the 
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addition of vectors in Z 3 , see BH93, Ch. 6] for more details. It is actually a graded K-algebra, with the 
grading being induced by the height n, i.e. by the decomposition C Pi Z 3 = (J n C n . If the variable Xk in A 
stands for the monomial s l P , we can identify it with the point 1) S C PI Z 3 . The multiplication of two 
monomials in A corresponds to the addition of two vectors in C . 

It is easy to verify that the above correspondence extends to a graded isomorphism of K-algebras between 
A and K[C] by observing that the relations of I(^) correspond to different decompositions of an element 
of C n as the sum of elements of smaller degree, so we actually have 

A ~ K[C] 

Exploiting this combinatorial description of the ring A we can deduce some algebraic properties, but let us 
first recall the definition of the canonical module, a notion we will use in this section. 

Definition 3. Let R = K[Xi, . . . , X n ], I an ideal of R and suppose that M = R/I is of dimension d. Then 
the canonical module of R is defined as lor = R[—n\ and 

lo m = Ext"~ d (M,R[-n\) 

is the canonical module u>m of M. 

The ring A is an affine n ormal semigroup ring by BH93, Prop. 6.1.2 and 6.1.4], since C H Z 3 is a normal 
semigroup. Moreover, by [BH93I . Prop. 6.3.5] it is Cohen-Macaulay and its canonical module u>a is the ideal 
generated by the monomials that correspond to integer points in the interior of C . This shows that A is 
Gorenstein if and only if the first Ci with non-empty interior (either i — 1, i = 2, or i — 3) contains exactly 
one point. In this case, it is actually easy to see the isomorphism between oja and A geometrically: It is 
nothing else th an the tr anslat ion th at mov es this point in the interior of Ci to the origin. Note that in the 
previous works [BJ03j . |BC05| . and |BD07j . the ring A was always Gorenstein and this property was used 



in some of the proofs. In our context we have to do without this property, which means that some of the 
proofs need to be modified. 

3. Homological tools 

3.1. Overview of approximation complexes 

W e will qu ic kly re call th e cons truction of the approximation complex Z, in order to fix notation, compare 
also |HSV83| . JVa9|, and |BJ03| . 



Let us denote by the class of the variable in the homogeneous coordinate ring A = K[2£]/J of 
where J = I{2f ) and X_ stands for the sequence Xo, ■ ■ ■ , X m . A is a graded ring, each variable having weight 
1. Let I = (.91,52,53,54) C A be the ideal generated by the gi, recall that d = deg(<?i). We consider the 
Koszul complex (K,(g, A), 5 m ) associated to 51, ... ,54 over A 

A[-Ad] A[-3d} A & A[~2d} 6 £ A[-d} 4 ^ A 

where the differentials are matrices with ±51,..., ±54 as non-zero entries. Write Kq = A, K\ = A^[— d], 
K 2 = A 6 [-2d], K 3 = A 4 [-3d], and K A = A[-id]. Set Z, := ker(^) C K i: which says that Z { also keeps the 
degree shift. Note that with this notation the sequence 

Zi -» Ki -» Bi..x -» (3) 

is an exact sequence of graded modules (with morphisms of degree zero). 

We set Zi = Zi[i ■ d] <£>aA[T^ : which we will consider as bigraded A[T]-modules (one grading is induced by 
the grading of A, the other one comes from setting deg(Ti) = 1 for all i). Now the approximation complex 
of cycles (Z,(g, A), e,), or simply Z», is the complex 

0^Z 3 (-3)^Z 2 (-2)^Z 1 (-l)^Z o 



(i 



where the differentials e. are obtained by replacing gi by Tj for all i in the matrices of 5, and where the 
degree shifts are with respect to the grading by the Ti. Then im(ei) is generated by the linear syzygies of 
the gi and 

H (Z.) = A[T]/im( ei ) ~ Sym A (I) 

From now on, when we take the degree v part of the approximation complex, denoted (-Z.)„, it should 
always be understood to be taken with respect to the grading of A. Hereafter we denote by m the maximal 
ideal (X ,...,X m ) CA. 

The geome tric in tuition behind the Z-complex is quite profound, we only give some hints and refer to ChOG, 



Sect. 3] or Va94| for a more thorough treatment of the subject. The symmetric algebra is closely related to 
the Rees algebra ReesA(-f), which can be defined as the quotient of A[T] by all the syzygies (not only the 
linear ones). One has thus a canonical surjection from Syrn A (I) onto Rees^/), which induces an inclusion 

Biproj(Rees A (/)) <-> Biproj(Sy mj4 (/)) (4) 

Now Biproj(ReesA(-0) corresponds to the closure of the graph of the map g and its image by the projection 
to P 3 equals the surface =5^, while Biproj(Sym j4 (/)) is a priori a bigger object. However, Sym A (I) is in 
some ways easier to study and under suitable conditions on the base points the inclusion in ((4J becomes 
an isomorphism and one can retrieve the information about 5? contained in the Rees algebra from the 
symmetric algebra. More precisely, we will see that the implicit equation of 5? can be obtained from the 
determinant of certain graded parts of the iJ-complex. 

The next lemma shows that the complex Z,(gi, . . . ,34; A) is acyclic if the base points are local complete 
inte rsectio ns and finite in number. This is a standard hypothesis for syzygy-based implicitization methods, 
KD06| . 
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Lemma 4. Let I = (31,32,33,34) C A. Suppose that 2? := Proj(A//) C has at most dimension and is 
locally a complete intersection, then the complex Z, is acyclic. 



Proof. If there ar e no b ase points this follows immediately from BJ03, Prop. 4.7] and for finitely many 



base points from [BJ03I . Prop. 4.9]. We only have to check that the hypotheses of these propositions are 
verified: In our case, we have n — 4 and we need to check that dim(A) = depth m (A) = n — 1 = 3, which 
is true because A is Cohen-Macaulay and because A is the homogeneous coordinate ring of a (projective) 
surface. Moreover, in the presence of base points, the equality depthj(A) = codim(I) = 2 = n — 2 is again 
a consequence of the Cohen-Macaulayness of A. □ 



Remark 5. It can be shown in a similar way as in [BD07I . Lemma 1] that the Z-complex is still acyclic if 



the base points are almost local complete intersections, but we will not treat this case here. 

3.2. Bounds on local cohomology 

The following lemma establishes a vanishing criterion on the local cohomology of Sym^ (/) , which ensures 
that the implicit equation c an be obtained as a generator of the annihilator of the symmetric algebra in a 



certain degree. We refer to [BS98I ] for more details on local cohomology, a detailed treatment of which is 
beyond the scope of this work. 

Lemma 6. Suppose that 3? := Proj(A//) C & has at most dimension and is locally a complete intersec- 
tion. If 7] is an integer such that 

H^(Sym A (I)) v = for all v > V 

then we have 

ann K[ T](Sym A (I) u ) = ann K[ T] (Sym A (I) v ) = kei(h) 

for all v > rj. 



Proof. The proof of BD071 Lemma 2] can be applied verbatim. □ 
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As we shall see, the annihilator in the above lemma can be computed as the determinant (or MacRae 
invariant) of the complex (Z m ) rp so we should give an explicit formula for the integer r\, but we first need 
to study the local cohomology of A using its combinatorial structure as a semigroup ring. The following 
definition is the same as [MS05I . Def. 11.15]. 

Definition 7. Let M be a graded ^-module. The Matlis dual M v of M is the ^-module defined by 

(M v )_ u =Hom K (M„,K), 

the multiplication being the transpose. One has (M v ) v = M if all the graded parts M u of M are finite- 
dimensional as K-vector spaces. 

Lemma 8. Let M be a finitely generated graded A-module of dimension r. Then M is Cohen- Macaulay if 
and only if H l m {M) = for all i ^ r and H^(M) = ufj^ is the Matlis dual to lom- 

Proof. This is [MS05I . Th. 13.37]. □ 

So the local cohomology of an A-module that is Cohen-Macaulay can be expressed in terms of its canonical 
module. Let us apply this to the A-module A. Using that dim(A) = 3 and that A is Cohen-Macaulay we 
immediately deduce 



Corollary 9. The local cohomology of A is 

H l m {A) 



if i ^ 3 
tv A if i = 3 



where u A is the Matlis dual to the canonical module uja- 

So the third local cohomology module of A is the only one that is non-zero. Actually, we do not need to 
know this module exactly; it is sufficient to know in which graded parts it vanishes. 

Corollary 10. Let a := max{i | C; contains no interior points} and let v£Z. Then we have H^A)^ = 
if v > ~a. 

Proof. By Corollary [9] and the definition of the Matlis dual we have the identities 

Hl{A) v = K)„ = Hom K ((w A )_ I/ ,K) 

but the module uja is generated by the elements in the interior of C, i.e. by elements of degree at least a + 1, 
so whenever v > ~a, it follows (uja)-v — and the modules in the above equation are all zero. □ 

We can now proceed to investigate the vanishi ng of th e Oth local cohom ology o f the symmetric algebra. The 
proof is similar to the corresponding theorems |BJ03l . 5.5 and 5.10] and |BD07 . Th. 1]. We give two bounds, 
an explicit one, which always holds, and a lower but more complicated bound for the case when there are 
base points. 

Theorem 11. Suppose that & := Proj(yl/7) C & has at most dimension and is locally a complete 
intersection. Then 

H^(Sym A (I)) v =0 Vv>v = 2d-a 

where a := max{i | C; contains no interior points } as before. Moreover, if there is at least one base point, 
one even has 

H^(Sym A (I)) v = > v Q = max{d - a, 2d + 1 - indeg(ff° (w a /Im> a ))}. 



Proof. The proof is virtually the same for the two cases. As the first one has been proven in Do08, Th. 
4.11], we only give a proof for the second bound. Consider the two spectral sequences associated to the 



double complex C^(Z t ), both converging to the hypercohomology of Z,. By Lemma SI -Z. is acyclic, hence 
the first spectral sequence stabilizes at step two with 



= ' 2 E% = K(H q (Z.)) 



fl*(Sym A (I)) for g = 0, 
otherwise. 



The second one has as first screen: 

'[El = Hl(Z q )[qd] ® A A[T lt . . .,T 4 ](-q). 

The comparison of the two spectral sequences shows that H^(Sym A (I)) v vanishes as soon as (i"E^) v 
vanishes for all p, in fact we have that 

end(H° m (Sym A (/))) < rqa«{end(i"£?)} = max{end«(Z p ) -p-d}. 

where we denote end(M) = maxjV | M„ / }. By Corollary [9] and the fact that Z$ = A, H^(Zq) = 0. 
Recall that the sequence 

-» -» -» Bi -f (5) 

shown in (|3|) is graded exact. From ([5]), applied to i = (writing £>o = 7) we obtain the long exact sequences 
of local cohomology 

Now H^(I) — 0, as I is an ideal of an integral domain, by Corollary O we have H^(Ki) = 0, hence H^(Zi) 
vanishes. 

Now as depth A (I) > 2, the Koszul complex is exact for i > 4 — 2 = 2, i.e. £^ = Z^. ft is clear by 
construction of the Koszul complex that Z 4 = and that B% = im^) ~ A[— d]. Using that H% l (A) v — for 
v > — a by Corollary [lOl we can deduce that H^(Z 3 ) U = H^{B^) V = if i> > d — a. It follows that 

for p = 0, f , or p > 3 
id! ,"£■;; i •■. ^ forp = 2 

f for p = 3 

It remains to determine e. From the short exact sequence — > £?i — > — > iJ; — > we get the exact sequence 

H^(Z 1 )^H^(H 1 )^H^(B 1 )^0, 

hence, as i?^,(Z 2 ) = H^(Bi) by ([5]), there is a surje ctive g raded map H^(Hi) -» iJ^(Z 2 ). 

Moreover, setting — * := Homgr A ( — ,A/m), by Ch04l . Lemma 5.8] we have the graded isomorphism 
(KiHt))* = K(H ( 9l , . . ,,g 4 ;u A ))[4d\ S ^//.ua)!^. Hence, we obtain 

endM) = end(tf£(Z 2 )[2c2]) 

<end(H^(H 1 ))-2d 

= -indeg(#£ (u a /Im> a )M) - 2d 

= 2d~mdeg{Hl(Lu A /I.Lu A )). 

We have shown that e < 2d — mdeg(H^(uj a / '1 'a)) , hence H^(Sym A (I)) v vanishes as soon as v > vq := 
max{d- a,2d+ 1 - indeg(7J° {ua/I-ua))}- □ 

Remark 12. Clearly, the advantage of the bound vq = 2d — a is that it does not require the computation 
of H^(lua/ I -loa)i which can turn out to be difficult even in simple examples. However, even thou gh it m ight 
not be obvious at first sight, the second bound is lower. For example, take the cas e studi ed in [BD07 |. i.e. 
N'(/) is a unit square and A is the quotient K[Xo, Xl, X 2 , X^I/XqX^ — X\Xi- By [BD07I . Prop. 2], we can 
identify lu a — A[-4 + 2], hence u Q — 2d+ 1 - indeg(i?£(A/I)[— 2]) = 2d- 1 - indeg(/ sat ), whereas the naive 




9 



bound would be 2d — a — 2d — 1. Simi larly, in the case = P 2 , our bound coincides with the known bound 
v = 2d - 2 - indeg(7 sat ) from [BC05I . Th. 3.2], as compared to 2d - 2. 

Also in the general case, one always has 2d+ 1 — indeg(_ff° (lua/I-^a)) <2d — a due to to a being generated 
in degree at least a + 1 as explained in Section |2~T1 and obviously d — a < 2d — a. 



4. The representation matrix 

It can now be deduced that the determinant of the ^.-complex is a power of the im plicit equation of ,5? . 
Indeed, using LemmaHJ LcmmaEJ and Theorem [TT| a completely analogous proof to BJ03L Th. 5.2] shows 
the following. 

Theorem 13. Suppose that & :— Proj(A//) C 3~ has at most dimension and is locally a complete 
intersection. Let a := max{i | C; contains no interior points} as before and Vq = 2d — a. For any integer 
v > vq the determinant D of the complex (Z») v of¥L[T]-modules defines (up to multiplication with a constant) 
the same non-zero element in K[T] and 

D = p de s(a) 

where F is the implicit equation of 5^ . 

By Theorem II 11 one can replace the bound in this result by the mor e preci se b ound vn = max{<i — a, 2d + 



ay iiieorem iin one can replace me uounu m uns result uy uic more precise uounu j/q — max^a — a, z,a -f 
1 — mdcg(H^(ujA/ Fuja))} if there is at least one base point. As in ChOtj ] or BCJ06| . there is a possible 



generalization of the above theorem to the case of almost local completion intersection base points. However, 



the proofs of the corresponding results (or the one of [BC05I . Th. 4]) do not apply directly here, because 
they use at some points that A is Gorenstein, which is not necessarily the case in the toric setting. 

By GKZ94I . Appendix A] , the determinant D can be computed either as an alternating sum of subdetermi- 



nants of the differentials in Z v or as the greatest common divisor of the maximal-size minors of the matrix 
M associated to the first map {Z\) v — * {Zq) v - Note that this matrix is nothing else than the matrix M v of 



linear syzygies as described in the introduction; it can be computed with the same algorithm as in [BD07 | 
by solving the linear system given by the degree vq part of ([1]). As an immediate corollary we deduce the 
following very simple translation of Thcorem ll31 which can be considered the main result of this paper. 



Corollary 14. Let 2? --■> P 3 be a parametrization of the surface y CP 3 given by g — (g\ : <?2 : 53 • 9i) 
with gi G A. Let M v be the matrix of linear syzygies of g%, . . . , 54 in degree v > 2d — a, i.e. the matrix of 
coefficients of a 'K-basis of Syz(g) l/ with respect to a WL-basis of A v . If g has only finitely many base points, 
which are local complete intersections, then M v is a representation matrix for the surface ,5f . 



We should also remark that by |KD06l . Prop. 1] (or [CoOll Appendix]) the de gree of the surface 5? can be 



expressed in terms of the area of the Newton polytope and the Hubert-Samuel multiplicities of the base 
points: 

deg( 5 )deg(^) = Area(N(/)) - £ e p (6) 

pev( gi ,..., g4 ,)cs- 

where Area(N(/)) is twice the Euclidean area of N(/), i.e. the normalized area of the polygon. For locally 
complete intersections, the multiplicity e p of the base point p is just the vector space dimension of the local 
quotient ring at p. 



5. The special case ST = P 1 x P 1 

Bihomogeneous parametrizations, i.e. the case J = P 1 x P 1 , are particularly important in practical ap- 
plications, so we will now make explicit the most important constructions in that case and make some 
refinements. We also include an implementation in Macaulay2 M2] in the Appendix. 
In this section, we consider a rational parametrization of a surface 5? 
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(s : u) x (t : v) i-> (/ x : f 2 : f 3 : U){s,u,t,v) 

where the polynomials fi, ■ . ■ , fi are bihomogeneous of bidegree (e±, e 2 ) with respect to the homogeneous 
variable pairs (s : u) and [t : v), and ei,e2 are positive integers. We make the same assumptions as in the 
general toric case. Let d — gcd(ei, e 2 ), e[ = % , and e' 2 = So we assume that the Newton polytope N(/) 
is a rectangle of length ei and width e 2 and N'(/) is a rectangle of length e[ and width e' 2 (in fact N(/) 
might be smaller, but in this section we homogenize with respect to the whole rectangle). 

So P 1 x P 1 can be embedded in P m , m = (e[ + l)(e' 2 + 1) — 1 through the Segre-Veronese embedding 

P = Pe lt e a 

P 1 x P 1 ^ P m 
(s : u) X (t : v) i-> (. . . : s'u^-'Pv^ : . . .) 

We denote by 3F its image, which is an irreducible surface in P m , whose ideal J is generated by quadratic 
binomials. We have the following commutative diagram. 

pi x pi_£^p3 (7) 
& 

with g — (gi : . . . : 34), the ^ being polynomials in the variables Xq, . . . ,X m of degree d. We denote by 
A = K[Xq, . . . , X m ]/J the homogeneous coordinate ring of & . We can give an alternative construction of 
the coordinate ring; consider the N-graded K-algebra 

S := ^ (K[s, u] ne > x ® K K[t, v} ne > 2 ) C K[s, u, t, v] 

which is finitely generated by 5i as an S'o-algebra. Then P 1 x P 1 is the bihomogeneous spectrum Biproj(S') 
of S, since Proj(0 neN K[s, u] ne ' ) = Proj(@ n6N K[t, v] ne > 2 ) = P 1 . The Segre-Veronese embedding p induces 
an isomorphism of N-graded K-algebras 

A^S 

where X l,J = X^ + iy i+ j for i = 0, . . . , e[ and j = 0, . . . , e' 2 and the implicit equation of .y can be obtained 
by the method of approximation complexes described in the previous sections as the kernel of the map 

K[T 1 ,...,T i ]^A 
Ti^g% 

The ring A is an affine normal semigroup ring and it is Cohen- Macaula y. It is Gorenstein if and only if 
e i = e 2 — 1 ( or equivalently e\ — e 2 ), which is the case treated in The ideal J is easier to describe 



than in the general toric case (compare |Su06l . 6.2] for the case e' 2 = 2). The generators of J can be described 
explicitly. Let 




X l ' e ' 2 

then the ideal J is generated by the 2-minors of the 4 x e^e^-matrix below built from the matrices Ai\ 

A ... A e ,^i 
A l ... A e , 



(8) 
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Let us also state the degree formula for this setting, which is a direct corollary of ©: 

deg(g)deg(^) = 2eie 2 e p 

pev( gi ,..., g4 )c3r 

where as before e p is the multiplicity of the base point p. 

We have claimed before that it is better to choose the toric variety defined by N'(/) instead of N(/). Let 
us now give some explanations why this is the case. As we have seen, a bihomogeneous parametrization of 
bidegree (ei, e^) gives rise to the toric variety 3" = P 1 x P 1 determined by a rectangle of length e[ and width 
e' 2 , where = ^f, d = gcd(ei, €2), and whose coordinate ring can be described as 

S := (K[s, u] ne; ® K K[t, v] ne , 2 ) c K[s, u, t, v] 

tiGN 

Instead of this embedding of P 1 x P 1 we could equally choose the embedding defined by N(/), i.e. a rectangle 
of length ei and width e2, in which case we obtain the following coordinate ring 

S := (J) (K[s,u] nei <Z>KK[t,v] ne2 ) C K[s,u,t,v] 
It is clear that this ring also defines P 1 x P 1 and we obviously have an isomorphism 

between the graded parts of the two rings, which means that the grading of S is coarser and contains 
less information. It is easy to check that the above isomorphism induces an isomorphism between the 
corresponding graded parts of the approximation complexes Z, corresponding to S and Z % corresponding 
to S, namely 

Z v — Zd. v 

If the optimal bound in Theorem [T^] for the complex Z is a multiple of d, i.e. i>q = d ■ 77, then the optimal 
bound for Z is i>q = rj and we obtain isomorphic complexes in these degrees and the matrix sizes will be 
equal in both cases. If not, the optimal bound Aj is the smallest integer bigger than ^ and in this case, 
the vector spaces in Zo will be of higher dimension than their counterparts in Z Vo and the matrices of the 
maps will be bigger. An example of this is given in the next section. 

6. Examples 

Example 15. We first treat some ex amples from [kD06]. Example 10 in the cited paper, which could not 
be solved in a satisfactory manner in [BD07], is a surface parametrized by 

fi = (t + t 2 )(s - l) 2 + (1 + st - s 2 t){t - l) 2 
h = i-t - t 2 )(s - if + (-1 + st + s 2 t){t - l) 2 
/ 3 = (t - t 2 )(s - l) 2 + (-1 - st + s 2 t)(t - l) 2 
U = (t + t 2 ){s - l) 2 + (-1 -st- s 2 t)(t - l) 2 
The Newton polytope N'(/) of this parametrization is 




1 2 
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We can compute the new parametrization over the associated variety, which is given by linear forms ffi, . . . , 54, 
i.e. d = 1 (since there is no smaller homothety N'(/) of N(/)) and the coordinate ring is A — K[Xq, . . . , Xg]/J 
where J is generated by 21 binomials of degrees 2 and 3. Recall that the 9 variables correspond to the 9 
points in the Newton polytope. In the optimal degree vq = 1 as in Theorem 1 1 1 1 the implicit equation of 
degree 5 of the surface 5^ is represented by a 9 x 14-matrix, compared to a 15 x 15-matrix with the toric 
resultant method (from which a 11 x 11-minor has to be computed) and a 5 x 5-matrix wit h the m ethod of 



moving planes and quadrics. Note also that this is a major improvement of the method in BD07J], where a 
36 x 42-matrix representation was computed for the same example. 

Example 16. Example 11 of |KD06j is similar to Example 10 but an additional term is added, which 
transforms the point (1,1) into a non-LCI base point. The parametrization is 

h = (t + t 2 )(s - if + (l + st- s 2 t)(t - if + (t + st + st 2 )(s - l)(t - 1) 
h = (-* - t 2 )(s - if + (-1 + st + s 2 t)(t - if + (t + st + st 2 )(s - l)(t - 1) 
f 3 = (i- t 2 )(s - if + (-1 -st + s 2 t)(t - if + (t + st + st 2 )(s - l)(t - 1) 
f 4 = (t + t 2 )(s - if + (-1 - st - s 2 t)(t - if + (t + st + st 2 )(s - l)(t - 1) 



The Newton polytope has not changed, so the embedding as a toric variety and the coordinate ring A are 
the same as in t he previous example. Again the new map is given by gi, ... ,54 of degree 1. 

As in KD06| . the method represents (with i/q = 1) the implicit equation of degree 5 times a linear 



extraneous factor caused by the non-LCI base point. While the Chow form method represents this polynomial 
as a 12 x 12-minor of a 15 x 15-matrix, our representation matrix is 9 x 13. Note that in this case, the method 
of moving lines and quadrics fails. 

Example 17. In this example, we will see that if the ring A is not Gorenstein, the correction term for i/q is 
different from indeg(/ sat ), unlike in the homogeneous and the unmixed bihomogeneous cases. Consider the 
parametrization 

h - (s 2 + t 2 )t 6 s 4 + (1 + s 3 t 4 - sH A ){t - i)V - 1) 
h = (s 2 - t 2 )t 6 s 4 + (-1 + s 3 t 4 + s 4 t 4 ){t - l) 5 (,s 2 - 1) 
h = (s 2 - t 2 )t 6 s 4 + (-1 - s 3 t 4 + s 4 t 4 )(t - 1) V - 1) 
/ 4 = (s 2 + t 2 )t 6 s 4 + (-1 - sh 4 - s 4 t 4 )(t - 1) V - 1) 

We will consider this as a bihomogeneous parametrization of bidegree (6,9), that is we will choose the 
embedding p corresponding to a rectangle of length 2 and width 3. The actual Newton polytope N(/) is 
smaller than the (6, 9)-rectangle, but does not allow a smaller homothety. One obtains A = K[Xo, . . . , Xu]/J, 
where J is generated by 43 quadratic binomials and the associated gi are of degree d = 3. It turns out that 
va = 4 is the lowest degree such that the implicit equation of degree 46 is represented as determinant of Z VQ , 
the matrix of the first map being of size 117 x 200. So we cannot compute vq as 2d — indeg(/ sat ) = 6 — 3 = 3, 
as one might have been tempted to conjecture based on the results of the homogeneous case. This is of 
course due to A not being Gorenstein, since the rectangle contains two interior points. 

Let us make a remark on the computation of the representation matrix. It turns out that this is highly 
efficient. Even if we choose the non-optimal bound v = 6 as given in Theorem 1131 the computation of the 
247 x 518 representation matrix is computed instantaneously in Macaulay2. Just to give an idea of what 
happens if we take higher degrees: For v = 30 a 5551 x 15566-matrix is computed in about 30 seconds, and 
for v = 50 we need slightly less than 5 minutes to compute a 15251 x 43946 matrix. 

In any case, the computation of the matrix is relatively cheap and the main interest in lowering the bound 
vq as much as possible is the reduction of the size of the matrix, not the time of its computation. This 
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reduction improves the performance of algorithmic applications of our approach, notably to decide whether 
a given point lies in the parametrized surface. 

Example 18. In the previous example, we did not fully exploit the structure of N(/) and chose a bigger 
polygon for the embedding. Here is an example where this is necessary to represent the implicit equation 
without extraneous factors. Take (/i, / 2 , / 3 , U) = (st 6 + 2, st 5 - 3st 3 , si 4 + 5s 2 i 6 , 2 + s 2 i 6 ). This is a very 
sparse parametrization and we have N(/) = N'(/). The coordinate ring is A = K.[Xo, . . . , X5]/ J, where 
,/ = (X 2 — X2X4, X2X3 — X1X4, Xf — X1X3, X\ — XqX^) and the new base-point-free parametrization g 
is given by (51, 32, 33; 9i) — (2X3 + X4, — 3Xx + X3, X2 + 5X 5 , 2X + X5). The Newton polytope looks as 
follows. 




For vq — Id = 2 we can compute the matrix of the first map of (Z m )^ , which is a 17 x 34- matrix. The 
greatest common divisor of the 17-minors of this matrix is the homogeneous implicit equation of the surface; 
it is of degree 6 in the variables T± , . . . , T4: 

2809T 1 2 T 2 4 + 124002T 2 6 - 5618T 1 3 T 2 2 T 3 + 66816TiT 2 4 T 3 + 2809T 1 4 T 2 
-50580T?t2t% + 86976T 2 4 T 3 2 + 212T 1 3 T 3 3 - 14210TiT 2 2 r 3 3 + 3078^^ 
+13632T 2 2 T 3 4 + 116TiT| + 841T 3 6 + 14045T 1 3 T 2 2 T 4 - 169849Ti T 2 4 T 4 
-14045T 1 4 T 3 T 4 + 261327T 1 2 r 2 2 T 3 T4 - 468288T 2 4 T 3 T 4 - 7208T 1 3 T 3 2 T 4 
+157155Tir 2 2 T 3 3 r 4 - 31098TfT|T 4 - 129215T 2 2 T 3 3 T 4 - 4528Ti T 3 4 T 4 
-12673T|T 4 - 16695T^T 2 2 T 4 2 + 169600T 2 4 T 4 2 + 30740rpT 3 T| 
-433384TiT 2 2 T 3 T 4 2 + 82434T 1 2 T 3 2 T 4 2 + 269745T 2 2 T 3 2 T 4 2 + 36696TiT 3 3 T 4 2 
+63946T 3 4 T| + 2775T x T|T| - 19470T 1 2 T 3 T| + 177675T 2 2 T 3 T 4 3 
-85360TiT 3 2 T 4 3 - 109490T 3 3 T 4 3 - 125T 2 2 T 4 4 + 2900Ti T 3 T 4 4 
+7325T 3 2 T 4 4 - 125T 3 T| 

As in Example [17] we could have considered the parametrization as a bihomogeneous map either of bidegree 
(2, 6) or of bidegree (1, 3), i.e. we could have chosen the corresponding rectangles instead of N(/). This leads 
to more complicated coordinate rings (20 resp. 7 variables and 160 resp. 15 generators of J) and to bigger 
matrices (of size 21 x 34 in both cases). Even more importantly, the parametrizations will have a non-LCI 
base point and the matrices do not represent the implicit equation but a multiple of it (of degree 9). Instead, 
if we consider the map as a homogeneous map of degree 8, the results are even worse: For vq = 6, the 
28 x 35-matrix M Vo represents a multiple of the implicit equation of degree 21. 
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To sum up, in this example the toric version of the method of approximation complexes works well, 



whereas it fails over 



and P~. This shows that the extension of the method to toric varieties really is 



a generalization and makes the method applicable to a larger class of parametrizations. 

Interestingly, we can even do better than with N(/) by choosing a smaller polytope. The philosophy is that 
the choice of the optimal polytope is a compromise between two criteria: 

- The polytope should be as simple as possible in order to avoid that the ring A becomes too complicated. 

- The polytope should respect the sparseness of the parametrization (i.e. be close to the Newton polytope) 
so that no base points appear which are not local complete intersections. 

So let us repeat the same example with another polytope Q, which is small enough to reduce the size of the 
matrix but which only adds well-behaved (i.e. local complete intersection) base points: 

3. 




^0 1 

The Newton polytope N(/) is contained in 2 • Q, so the parametrization will factor through the toric variety 
associated to Q, more precisely we obtain a new parametrization defined by 

{9U92, 33, 34) = (2X 2 + SXoXi + X 2 X 4 , X X X 4 + 5Xl 2X 2 + X\) 

over the coordinate ring A = K[X 0) ■ ■ . ,X 4 ]/J with J = (Xf - X 1 X 3 ,X 1 X 2 - X X 3 ,Xf - X X 2 ). The 
optimal bound is Vq — 2 and in this degree the implicit equation is represented directly without extraneous 
factors by a 12 x 19-matrix, which is smaller than the 17 x 34 we had before. 

Example 19. As we have seen, the size of the matrix representation depends on the given parametrization 
and as a preconditioning step it is often advantageous to cho ose a simpler parametrization of the same surface, 
if that is possible. For example, approaches such as [Sc03j | can be used to find a simpler reparametrization 
of the given surface and optimize the presented methods. 

Another important factor to consider is that all the methods we have seen represent the implicit equation 
to the power of the degree of the parametrization. On one hand, it can be seen as an advantage that this 
piece of geometric information is encoded in the matrix representation, but on the other hand, for certain 
applications one might be willing to sacrifice the information about the parametric degree in order to obtain 
smaller matrices. If this is the case, there exist (f or cert ain surface parametrizations) algorithms to compute 
a proper reparametrization of the surface, e.g. |Pe06| . and in these cases it is highly advisable to do so 
before computing the matrix representation, because this will allow us to represent the implicit equation 
directly instead o f one of its powers, and the matrices will be significantly smaller. Let us illustrate this with 
Example 2 from *Pe06| , which treats a parametrization / defined by 

fi = {sH A + 2s A t 2 + 5s 4 + 2t A + 4t 2 + ll)(.s 4 + 1) 
h = (s 4 t 4 + 2sH 2 + 5s 4 + t 4 + 2t 2 + 6) 
h = -(s 4 t 4 + 2sH 2 + 5.s 4 + t 4 + 2t 2 + 3)(s 4 + 1) 
/ 4 = (i 4 + 2< 2 + 5)( S 4 + l) 

This is a parametrization of bidegree (8, 4) and its Newton polytope is the whole rectangle of length 8 and 
width 4, so we can apply the method of approximation complexes for P 1 x P 1 . We obtain a matrix of size 
45 x 59 representing Fij?, where 



= 2T{T 2 - T 2 T 3 - 3T1T4 - 2T 2 T 4 + 3Tf 
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is the implicit equation and deg(/) = 16. Using the algorithm presented in Pc06] one can compute the 
following proper reparametrization of the surface 5?: 



/i = -(11 + st - 5s - 2t)(s - 1) 
/ 2 = 6 - 1 - 5s + st 
f 3 = (-t + st-5s + 3)(s-l) 
A = (t-5)(a-l) 

This parametrization of bidegree (2, 1) represents Fy directly by a 6 x 11-matrix. 
7. Final remarks 

Representation matrices can be efficiently constructed by solving a linear system of relatively small size (in 
our case dim^A^+d) equations in Adim^(A^) variables). This means that their computation is much faster 
than the computation of the implicit equation and they are thus an interesting alternative as an implicit 
representation of the surface. 

In this paper, we have extended the method of matrix representations by linear syzygies to the case of ra- 
tional surfaces parametrized over toric varieties (and in particular to bihomogeneous parametrizations) . This 
generalization provides a better understanding of the method through the use of combinatorial commutative 
algebra. From a practical point of view, it is also a major improvement, as it makes the method applicable 
for a much wider range of parametrizations (for example, by avoiding unnecessary base points with bad 
properties) and leads to significantly smaller representation matrices. Let us sum up the advantages and 
disadva ntages compared to other techniques to compute matrix representations (e.g. the ones introduced in 
The most important advantages are: 

- The method works in a very general setting and makes only minimal assumptions on the parametrization. 
In particular, it works well in the presence of base points. 

- Unlike the method of toric resultants, we do not have to extract a maximal minor of unknown size, since 
the matrices are generically of full rank. 

- The structure of the Newton polytope of the parametrization is exploited, so one obtains much better 
results for sparse parametrizations, both in terms of computation time and in terms of the size of the 
representation matrix. Moreover, it subsumes the known method of approximation complexes in the case 
of dense homogeneous parametrizations, in which case the methods coincide. 

Disadvantages of the method are the following. 

- Unlike with the toric resultant or the method of moving planes and surfaces, the matrix representations 
are not square. 

- The matrices involved are generally bigger than with the method of moving planes and surfaces. 

It is important to remark that those disadvantages are inherent to the choice of the method: A square matrix 
built from linear syzygies does not exist in general and it is an automatic consequence that if one only uses 
linear syzygies to construct the matrix, it has to be bigger than a matrix which also uses entries of higher 
degree. The choice of the method to use depends very much on the given parametrization and on what one 
needs to do with the matrix representation. 



Appendix: Implementation in Macaulay2 In this appendix we show how to compute a matrix represen- 
tation with the method developed in this paper, using the computer algebra system Macaulay2 [M2j . As 
it is probably the most interesting case from a practical point of view, we restrict our computations to 
bi-homogeneous parametrizations of a certain bi-degree (e±, e%). However, the method is easily adaptable to 
the toric case, or more precisely to a given fixed Newton polytope N(/) and, where it is appropriate, we 
will give hints on what to change in the code. Moreover, we are not claiming that our implementation is 
optimized for efficiency; anyone trying to implement the method to solve computationally involved examples 
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is well-advised to give more ample consideration to this issue. For example, in th e toric case there are better 
suited software systems to compute the generators of the toric ideal J, see 4ti2| . 

Let us start by defining the parametrization / given by (/i, . . . , f^). 
S=QQ[s ) u ) t,v] ; 
el=4; 
e2=2; 

f l=s~4*t~2+2*s*u~3*v~2 
f 2=s~2*u~2*t*v-3*u~4*t*v 
f 3=s*u~3*t*v+5*s~4*t~2 
f 4=2*s*u~3*v~2+s~2*u~2*t*v 
F=matrix{{f 1 , f 2 , f 3 , f 4}} 

The reader can experiment with the implementation simply by changing the definition of the polynomials 

and their degrees, the rest of the code being identical. We first set up the list st of monomials s % V of bidegree 

(e'i, e' 2 ). In the toric case, this list should only contain the monomials corresponding to points in the Newton 

polytope N'(/). 

st={}; 

1=-1; 

d=gcd(el,e2) 
eel=mimerator (el/d) ; 
ee2=numerator (e2/d) ; 

for i from to eel do ( 
for j from to ee2 do ( 

st=append(st , s~i*u~ (eel-i) *t~ j*v~ (ee2-j ) ) ; 
1=1+1 

) 

) 

We compute the ideal J and the quotient ring A. This is done by a Grobner basis computation which works 
well for examples of small degree, but which should be replaced by the matrix formula_in ([8]) for more 
complicated examples. In the toric case, there exist specialized software systems such as [4ti2| to compute 
the ideal J. 

SX=QQ [s , u , t , v , w , x_0 . . x_l , MonomialOrder=>Eliminate 5] 
X={}; 

st=matrix {st}; 

F=sub(F,SX) 

st=sub(st,SX) 

te=l; 

for i from to 1 do ( te=te*x_i ) 

J=ideal (l-w*te) 

for i from to 1 do ( 

J=J+ideal (x_i - st_(0,i)) 

) 

J= selectlnSubringd ,gens gb J) 

R=QQ[x_0. .x_l] 

J=sub(J,R) 

A=R/ideal(J) 

Next, we set up the list ST of monomials s l P of bidegree (e\, e<i) and the list X of the corresponding elements 
of the quotient ring A. In the toric case, this list should only contain the monomials corresponding to points 
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in the Newton polytope N(/). 

use SX 

ST={>; 

for i from to el do ( 

for j from to e2 do ( 

ST=append(ST, s~i*u~ (el-i) *t~ j*v~ (e2-j ) ) ; 
) 

) 

X={}; 

for z from to length(ST) -1 do ( 
f=ST_z; 
xx=l ; 

is=degree substitute (f , -[u=>l , v=>l ,t=>l}) ; 
is=is_0 ; 

it=degree substitute (f , -[u=>l , v=>l , s=>l}) ; 
it=it_0; 

iu=degree substitute (f , {t=>l , v=>l , s=>l}-) ; 
iu=iu_0; 

iv=degree substitute (f , {u=>l ,t=>l , s=>l}-) ; 
iv=iv_0; 
ded=0; 
while ded < k do ( 

for mm from to 1 do ( 

js=degree substitute (st_ (0 , mm) , {u=>l , v=>l ,t=>l}) ; 
js=js_0; 

jt=degree substitute (st_ (0 ,mm) , {u=>l , v=>l , s=>l}) ; 
jt=jt_0; 

ju=degree substitute (st_ (0 ,mm) , {t=>l , v=>l , s=>l}) ; 
ju=ju_0; 

j v=degree substitute (st_ (0 ,mm) , {u=>l ,t=>l , s=>l}) ; 
jv=jv_0; 

if is>=js and it>=jt and iu>=ju and iv>=jv then ( 

xx=xx*x_mm; 

ded=ded+l; 

is=is-j s ; 

it=it-jt ; 

iv=iv-jv; 

iu=iu-ju; ))); 
X=append(X,xx) ; ) 
We can now define the new parametrization g by the polynomials g\, ... ,54. 
X=matrix {X} ; 
X=sub(X,SX) 

(M , C) =coef f icients (F , Variables=> 

{s_SX , u_SX , t_SX , v_SX> , Monomials=>ST) 

G=X*C 

G=matrix{{G_(0,0) ,G_(0,1) ,G_(0,2) ,G_(0,3)}} 
G=sub(G,A) 

In the following, we construct the matrix representation M. For simplicity, we compute the whole module 
Z\, which is not necessary as we only need the graded part [Z\) Vo . In complicated examples, one should 
compute only this graded part by directly solving the linear system given by ([TJ in degree vq. Remark that 
the best bound nu — v$ depends on the parametrization. 
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use A 

Zl=kernel koszul(l,G); 

nu=2*d-l 

S=A [T1,T2,T3,T4] 

G=sub(G,S) ; 

Zlnu=super basis (mi+d, Zl) ; 

Tnu=matrix{{Tl,T2,T3,T4}}*substitute(Zlnu,S) ; 

lll=matrix {{x_0..x_l» 

lll=sub(lll,S) 

11=0 

for i from to 1 do { ll=append(ll,lll_(0,i)) > 
(m,M)=coef f icients(Tnu,Variables=> 

ll,Monomials=>substitute(basis(nu,A) ,S)) ; 

M; 

The matrix M is the desired matrix representation of the surface 5? . 



Acknowledgements We thank Laurent Buse and Marc Chardin for useful discussions. 
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